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D ■ Abstract: Some supersymmetric models predicted long lived charged 

Q ■ sleptons. Following Albuquerque, Burdman and Chackopp, we consider 

■^j- | an exact process where the r, the supersymmetry partner of r lepton, are 

produced inside the earth by collisions of high energy neutrinos with nucleons. 
We detailedly investigate the possible signals of upgoing r by comparing with 
C<") ! the background muons where atmospheric neutrino flux is taken into account. 

' The realistic spectra shows that km-scale experiments could see as many as 

■ 69 events a year by using the Waxman-Bahcall limit on the extraterrestrial 

neutrino flux. 
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1 Introduction 



Substantial evidences exist suggesting that most of the Universe's matter is non 
^ | luminous [2] EH- There are many predictions for its composition but the nature of dark 

$_i ' matter is yet unknown. 

c3 ; 

Weak scale supersymmetric theories of physics beyond the standard model, provide 
perhaps the most promising candidates for dark matter[3]. However, supersymmetry 
must be broken at the energy scale accelerator having reached since the super part- 
ners have not been observed yet. Supersymmetric models typically have a symmetry, 
called R-parity, which exclusively ensures that the Lightest Supersymmetric Particle 
(LSP) is the most stable. Obviously, the LSP is the natural candidate for dark matter. 
Different scales of supersymmetry breaking determine what the LSP is. Typically, if 
supersymmetry is broken at high scales, the LSP is likely to be the neutralino or the 
quintessino[3], however, if supersymmetry is broken at lower scales, the LSP is likely 
to be the gravitinojH]. In the models where the LSP is typically the quintessino or the 
gravitino, the Next to Lightest Supersymmetric Particle (NLSP) tends to be a charged 
slepton, typically the right-handed r, which has a long lifetime between microsecond 
and a year aroundj^j 0; depending on the scale of supersymmetry breaking and the 
slepton's mass. Of course, these lifetimes are negligible comparing with the age of our 
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universe, therefore, almost all these NLSP produced in the evolutive history of the 
universe decayed into LSP. 

Lest to mislead the reader, we should clarify several points. It is not absolutely cer- 
tain that the NLSP is either charged slepton or other neutral supersymmetric particle; 
furthermore, there is at present no direct evidence for the existence of supersymmetry. 
These are still unproven ideas. But some theoretical conclusions prefer the slepton to 
other neutral NLSPs [ZJ- Although speculative, supersymmetric dark matter is very 
well motivated and based on a simple physical principle. 

In the models of Slepton as NLSP, since the slepton is massive, actually, it can 
possibly be produced by high energy process. Collisions of high energy neutrinos with 
nucleons in the earth at energies above threshold for supersymmetric production can 
produce supersymmetry pairs which are unstable and promptly decay into slepton; 
since the slepton is charged and long lived, it is worthy to see whether its upgoing 
tracks could be detected by some large cosmic ray experiments or neutrino telescopes, 
such as L3+C (Sj, ICECUBE [9 a , Super-K |10j . etc. Based on a gravitino-LSP scenario, 
Albuquerque et al initiated that one can take neutrino telescopes as a direct probe of 
supersymmetry breaking^. Motivated by their work, our following calculations base 
on the quintessino-LSP scenario from Ref.jSj where it restricts the stau mass between 
WOGeV and ITeV and the lifetime between 10 6 ~ 10 7 seconds. 

2 Fundamental Process Analysis 

High energy neutrinos interacting with nucleons will produce supersymmetric par- 
ticles, and that will promptly decay into a pair of NLSPs, which have long life- 
time. This lifetime is definitely large enough so that r do not decay inside the earth. 
vN — > Iq —¥ 2r, the dominant process is analogous to the standard model charged cur- 
rent interactions, i.e, v^N — > u~ + anything. We include the corresponding processes 
at the parton level in Figure ^ 

We make use of the stau mass = 143GeV and two values for the squark masses: 
nriq = 150, 300GeV jllj in our calculations, and the cross section for supersymmetric 
production as a function of the neutrino energy is given in Figure |2j Also plotted for 
comparison is the standard model charged current cross section ^2] • From Figure |2] 
one can conclude supersymmetric interactions is much weaker than that of standard 
model. 

3 Calculate upgoing fluxes 

In order to be comparable with the results form Ref.pQ, we consider the similar 
physical process. We take the earth as the target and consider the incoming of some 
diffuse fluxes of high energy neutrinos. In our calculations, we make use of a model of 
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Figure 1: Feynman diagrams for supersymmetric particles production in uN collisions, 
the upper two diagrams refer to charged current interactions, and the under two dia- 
grams account for neutral current process, with the exchanges of the chargino x + an d 
neutralino x°, respectively. 
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Figure 2: The cross section of vN interactions as a function of neutrino energy. The 
curves correspond to = 143GeV and for squark masses rriq = 150,300GeF. The 
top curve refers to the standard model charged current interactions. 
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the earth density profile as detailed in Ref.|12j. 

Dynamic analysis shows that the energy threshold for r production has to be over 
~ lOOTeV. For atmospheric neutrino flux, usually experimental data give its ranges 
from 1 to lOTeV, which seems out of our need, however, theoretical results show us 
that atmospheric neutrino spectrum with energy above WTeV still act an important 
flux for t and muon production. We make use of atmospheric neutrino flux with energy 
under WTeV in Ref.|14j. and the flux with energy above WTeV in Ref. |13j. 

Figure 01 summarizes some popular theoretical and experimental upper limits on 
diffuse neutrino fluxes. The Waxman-Bahcall (WB) upper bound |15| assumes that 
100% of the energy of cosmic ray protons are lost to tt + and ir~ and that the tt + 
all decay to muons that also produce neutrinos. Ref. [To] also discussed the maximum 
contribution due to possible extra-galactic component of lower-energy < 10 17 eV, where 
protons have been first considered (max. extra-galactic p). Experimentally, AMANDA 
experiment gave a upper bound on diffuse neutrino flux |16j . By considering optically 
thick AGN models or by involving very strong magnetic fields, Mannheim, Protheore, 
and Rachen (MPR) have argued that one might be able to avoid the WB limit and get a 
higher upper bound JJj . From Figure |H1 actually, for extraterrestrial neutrino flux, we 
can conclude that the upper bound has been decided by AMANDA experiment on the 
neutrino energy below 10 e GeV , therefore, we take "max. extra-galactic p" as supplement 
for the neutrino energy above W^GeV. In this article, unless extra specification, the 
following results all make use of atmospheric neutrion flux with energy above lGeV 
adding the conservative extraterrestrial WB Limit as the incoming neutrino flux. 

There are two factors must be taken into account: one is the attenuation of high 
energy neutrinos through the earth, the other is the energy loss of charged r before 
running into a detector. 

The earth is opaque to ultra-high energy neutrinos. When the energy of neutrinos 
is over 40TeV, the vN interaction length turns out to be larger than the diameter of 
the earth, the attenuation of neutrinos must be considered. 

l„(E u ,x)=l°(E u )eM-x/l), (1) 

where I®(E U ) is the initial incoming neutrino flux, x is the depth a neutrino penetrat- 
ing the earth, and I refers to the interaction length, / = \j{otot ■ n), n accounts for the 
number density of the medium. Here crtoti i n principle, include all charged and neutral 
current processes contributing by both standard model and supersymmetric produc- 
tion. Since the initial interactions produce slepton and these are nearly degenerated in 
flavor, the flavor of the initial neutrino does not affect the results. 

The upgoing r event rate depends on the vN cross section in two ways: through the 
interaction length which governs the attenuation of the neutrino flux due to interactions 
in the earth, and through the probability that the neutrino converts to a r energetic 
enough to arrive at the detector with large than the threshold energy E™ m . The 
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Figure 3: Neutrino fluxes given by different models. The WB bound line gives the 
upper bound corrected for neutrino energy loss duo to redshift and for the maximum 
known redshift evolution. The dot curve is the maximum contribution due to possible 
extra-galactic component of lower-energy. The dash-dot line shows the experimental 
upper bound on diffuse neutrino flux established by the AMANDA experiment. The 
top line is given by authors Mannheim, Protheore, and Rachen (MPR) with considering 
optically thick AGN models or very strong magnetic fields. Dashed line is the theoretic 
predict for energetic atmospheric neutrino flux. 
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probability that a r produced in a vN interaction arrives in a detector with an energy 
above the r energy shreshold E~ un depends on the range R of a r in rock, which follows 
from the energy-loss relation |18j 

-&Er/dx = a + Ef/£, (2) 

here, the coefficients a and £ characterize the ionization and radiation losses respec- 
tively. For numerical estimates of ionization loss here we use a = 2MeV/(gcm- 2 ) [TO] . 
For a given momentum impulse, the radiation energy loss is inversely proportional to 
the square of the mass of the radiation particle. Thus the radiation length for r is ap- 
proximately {m^/rrifj) 2 times large than that for muons. We take £^ ~ 2.5 x 10 5 <?/cm 2 
in rock ^Bl and define e = a£, then the range can be expressed 

R(E~ T ,Er n )=^^^-. (3) 

The general solution of equation © is 

= (e'~ + e) exp(-x/0 - e. (4) 

The left side of (@J is to be interpreted as the residual energy of a r of initial energy 
E~ after penetrating a depth x of material. 

For those r ranging into the detector and fixing energy at E^, they can be con- 
tributed by any initial r with the energy above E^ and being produced at the distance 
R(E~, E^), therefore, the differential flux intensity can be expressed: 

dN~ l' R I' d(T~ 

—L=2tt sinOde / N A p(r)dx l v {E v1 x') — r dE u , (5) 

dE 7 J Jo J E th dE u dE~ 

Here, we define the zenith angle 6 as the angle between the incident direction of neu- 
trinos and the direction of the line linking the center of both the earth and the de- 
tector. Na is Avogadro's number, p(r) corresponds to the earth density, and r = 



RL + 2xRq) cos 0) is the distance from the center of the earth, where Rq refers 
to the radius of the earth. E^ is neutrino threshold energy for r productions. — d<J? , 

refers to the differential cross section of r productions. T u (E u ,x') and E~ are given by 
Equation Q and @, respectively. We take x' = 2R^ cos 8 — x instead of the distance 
a neutrino travelling in the earth. 



4 Possible signals analysis 

High energy vN processes can produce r as well as muons. Obviously, the upgoing 
muons, as the background flux, will range into the detector accompanying with r by 



6 




Figure 4: Energy spectrum of r pair events per km 2 , per year, for nig = 150, 300GeV. 
Also shown are the upgoing neutrino and muon flux through the detector. We make 
use of atmospheric neutrion flux with energy above lGeV adding the conservative 
extraterrestrial WB Limit as the incoming neutrino flux. 
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Figure 5: Angular distribution of f pair events per km 2 , per year, for rriq = IhQGeV . 
Also shown are the muon flux. We make use of atmospheric neutrion flux with en- 
ergy above lGeV adding the conservative extraterrestrial WB Limit as the incoming 
neutrino flux. 

similar interacting process. There are several possible ways exist to distinguish the r 
signals from the background muons. 

The first is its differential energy spectrum. Some detectors have very well energy 
resolution, typically, L3+C detector [Hj, can reconstruct the exact tracks of charged par- 
ticles in the magnetic field, which can determine the momentum, charges and direction 
of the incident particles at last. This can possibly provides the direct evidences for 
some exotic particles. In Figure^jwe show the energy distribution for the r pair events 
for two choices of q masses: 150GeV, 300GeV . Also shown is the upgoing neutrino 
flux as well as the energy distribution of upgoing muons'. We see that, the dominant 
contribution of r comes from the energy zones between 10 5 ~ W 7 GeV. This mainly 
because most of f are produced in the earth with the energy above 10 5 GeV and range 
into detector with high energy as well. However, with the energy of neutrinos growing, 
the range of r with the energy above about W 8 GeV turns out to be lager than the 
diameter of the earth, therefore, all initial r can totally arrive at the detector, which 
makes the events rate curve almost parallel with that of incoming neutrinos. Comparing 
with the curve of r flux, we can see, the muon flux's is much stronger at the low en- 
ergy zones(< 10 5 GeV), and turn out to be weaker at the high energy zones(> 10 5 GeV). 

The second is the angular distribution of events rate. Generally, neutrino telescopes 
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Figure 6: The Monte Carlo result for integral total 1000 pair events vs. their distance 
5L. We make use of atmospheric neutrion flux with energy above lGeV adding the 
conservative extraterrestrial WB Limit as the incoming neutrino flux. 

have poor energy resolution, only the angular resolution is used to reduce the back- 
ground. Figure [3 shows us the angular distribution of events rate of r and muons. From 
the figure one can conclude that the angular distribution of upgoing r is definitely differ- 
ent from muon's. The dominant contribution of upgoning r comes from the solid angle 
between 70 ~ 80 degrees. This is mainly because most of r have their energy between 
10 5 ~ W 7 GeV which hold the range about 10 cmwe, this approximately equal to the 
acclivitous thickness of the earth at the directions between 70 ~ 80 degrees. Therefore, 
to the isotropic incoming neutrinos within that directions, once interact with nucleons, 
the produced r will range into the detector totally. For most of muons, their ranges 
are neglectable comparing with the diameter of the earth, so the events will grow as 
the angle increasing. It is a pity that one can hardly reduce background through this 
method duo to the peak signals of angular distribution completely overlayed by muons'. 
Comparing with extraterrestrial neutrinos, atmospheric neutrino flux gives dominant 
contribution to the muon production. 

The third is the tracks of the r. The NLSPs are produced in pairs and that will 
promptly decay into r with the average energy {E~)i, 2 ~ 70%(20%)£„ [LI]. As men- 
tioned in Ref.pP, typical signal events include two tracks separated by certain distance 

5L m D&, (6) 
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Table 1: Number of events per km 2 per year for taking WB, AMANDA+max. extra- 
galactic p and MPR as extraterrestrial fluxes, respectively. Atmospheric neutrino flux 
with the energy above lGeV is taken into account as well. The first column refers to 
upgoing muons. The last two columns correspond to upgoing r for two different choices 
of squark masses: 150, 300Gey. 





muon 


m 4 = IbOGeV 


rriq = 300GeV 


Atmos.+ WB 


4.83 x 10 5 


69 


26 


Atmos.+ AMANDA+max. 


4.90 x 10 5 


470 


139 


Atmos.+ MPR 


5.02 x 10 s 


3164 


1244 



where D refers to the distance to the production point, and $ is the angle between 
the pair particles. Considering D to be the same order with r range, typically several 
hundred kilometers and i? within 10 -3 , 5L should typically be within several hundred 
meters. Figure El gives a Monte Carlo result for 1000 r pair events. We make use of at- 
mospheric neutino flux adding extraterrestrial neutrino flux(WB) as input and consider 
an exact interacting process inside the earth, including different incoming directions, 
attenuation of neutrinos and r, enery conversion from a neutrino to two r, the range 
of a t running in the earth, etc. At last we get an integral distribution for number 
of pair events(< 5L, Wm/bin) vs. the distance 5L, from which one can conclude that 
typical pair events hold the 5L about 120m(~ 50%) and most of them are within the 
confines of 320m(~ 80%). Owing to all muon signals are single tracks, seeking double 
tracks signals turns out to be an effective method to distinguish from background, if 
we take the two tracks which both enter into a detector within 6t = 4 microsecond 
as a double-track event, then for per km 2 per year, the accidental coincidence rate of 
muons is as few as 2N^N^5t ~ 6 x 10~ 2 , where refers to the number of events for 
muons. 

In Tabled we show the events rate for r pair production per km 2 per year on the 
atmospheric neutrino flux with the energy above lGeV adding different extraterrestrial 
neutrino fluxes. For being comparable, we also show the rates of upgoing muons. 
Comparing with the results from Ref . , the number of events are increased at least an 
order of magnitude. This mainly because we work under the SUSY breaking scenario 
with a quintessino LSP instead of a gravitino LSP. According to the analysis given 
above, some large scaled detectors appear to be sensitive to the relatively long lived r. 

5 Conclusions 

We discussed the exact production process of upgoing sleptons based on the long 
lifetime scenario with a quintessino LSP. Furthermore, we detailedly investigated the 
possible signals of events: energy spectrum, angular distribution and their tracks, from 
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which one can conclude that the characteristic signals of events are distinctively dif- 
ferent from that of background muons'. The study shows that seeking double-track 
upgoing events is a feasible method to detect r signals. It's worth to mention that 
we introduced the contribution from atmospheric neutrino flux and found it's very 
important for upgoing r and muon production. The event rates are also given, the 
quintessino-LSP scenario predicts more NLSPs than that of gravitino-LSP. The nu- 
meric results show that km-scale detectors are hopeful to get some positive results. 

The possible signals may provide a direct evidence for supersymmetry theory, fur- 
thermore, it can also offer a potential solution for dark matter problem. 
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